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Abstract The exact analytic solution is introduced for the rotational motion 
of a rigid body having three equal principal moments of inertia and subjected 
to an external torque vector which is constant for an observer fixed with the 
body, and to arbitrary initial angular velocity. In the paper a parametrization 
of the rotation by three complex numbers is used. In particular, the rows of 
the rotation matrix are seen as elements of the unit sphere and projected, by 
stereographic projection, onto points on the complex plane. In this represen- 
tation, the kinematic differential equation reduces to an equation of Riccati 
type, which is solved through appropriate choices of substitutions, thereby 
yielding an analytic solution in terms of confluent hypergeometric functions. 
The rotation matrix is recovered from the three complex rotation variables 
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by inverse stereographic map. The results of a numerical experiment con- 
firming the exactness of the analytic solution are reported. The newly found 
analytic solution is valid for any motion time length and rotation amplitude. 
The present paper adds a further element to the small set of special cases for 
which an exact solution of the rotational motion of a rigid body exists. 

Keywords rigid body dynamics, kinematics, rotation, integrable cases of 
motion, spherical ellipsoid of inertia 

1 Introduction 

The rotational motion of a rigid body has been studied for centuries. Besides 
the theoretical interest, the subject is of practical importance for many ap- 
plications, pertaining, in particular, to astronautics and celestial mechanics. 

The problem of the rotational motion of a rigid body can be divided into 
two parts. The dynamic problem aims to obtain the angular velocity of the 
rigid body with respect to an inertial reference by starting from the knowl- 
edge of the initial angular velocity and the history of the applied torque. On 
the other hand, the kinematic problem focuses on determining the current 
orientation of the body from the knowledge of the initial orientation and the 
history of the angular velocity. The differential equations governing the rota- 
tional dynamic and kinematic problems are decoupled if the applied torque 
does not depend on the body orientation. 

The exact analytic solution for the rotational motion of a rigid body exists 
only for a small number of special cases. 

When no torque acts on the body, exact solutions exist for both the 
dynamic and kinematic problems in the case of asymmetric rigid body with 
a fixed point coincident with the center of mass (Euler-Poinsot case) , and in 
the one of axis-symmetric rigid body under gravity force with a fixed point 
along the axis of symmetry (Lagrange-Poisson case) (Golubev 1960; Leimanis 
1965; Morton et al. 1974). 
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On the contrary, when the rigid body is subjected to external torques, the 
exact analytic solution exists, limited to the the dynamic problem, only for 
an axis-symmetric rigid body under constant axial torque (Bodcwadt 1952; 
Longuski 1991). Nevertheless, no exact analytic solution has been known for 
the kinematic problem, with the exception of the simple case of torque along 
a principal axis of inertia and initial angular velocity along the same axis. 

Several researchers have proposed approximate solutions. In particular, 
Longuski 1991, and Tsiotras and Longuski 1991, and Longuski and Tsiotras 
1995 use asymptotic and series expansion to get the approximate solution 
of the dynamics and kinematics for both the constant and the time-varying 
torque cases. Livneh and Wie 1997 analyze qualitatively, in the phase space, 
the dynamics of a rigid-body subjected to constant torque. As regards the 
kinematic problem, Iserles and N0rsett 1999 study the solution in terms of 
series expansion for the more general problem of solving linear differential 
equation in Lie Groups, building upon the work of Magnus 1954. Finally, 
Celledoni and Saefstroem 2006 propose ad-hoc numerical integration algo- 
rithms. 

For the first time, to the knowledge of the author, the present paper in- 
troduces the exact analytic solution for the rotational kinematic problem of 
a rigid body with a spherical ellipsoid of inertia subjected to a torque vector 
which is constant with respect to an observer fixed with the body, and to 
arbitrary initial angular velocity Besides a uniform sphere, any uniform reg- 
ular simplex has spherical ellipsoid of inertia (Klamkin, 1993). The proposed 
analytic solution is valid for any length of time and rotation amplitude, as 
opposed to the previously known approximate solutions cited above. 

In particular, the parametrization of the rotation matrix in terms of 
complex stereographic rotation variables, introduced earlier by Tsiotras and 
Longuski 1991, is used in the paper in order to obtain a kinematic differential 
equation of Riccati type which is then solved analytically by substitutions of 
variables. The solution involves the confluent hypergeometric function. 
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The proposed analytic solution is of considerable theoretical interest, as 
it adds a new element to the set of few special cases for which an exact 
solution of the rotational motion of a rigid body exists. Therefore, it can 
be used to rapidly obtain accurate solutions over arbitrary large intervals of 
time, and constitutes a significant new comparison case for error analysis of 
approximate algorithms. 

The paper is organized as follows: Section 2 introduces the dynamic and 
kinematic equations for the rotational motion of a generic rigid body, and 
the parametrization of the rotation by stereographic projection. Section 3 
reports the equations for the particular case of a rigid body having spherical 
ellipsoid of inertia and subjected to a constant torque. Section 4 introduces 
the new analytic integration of those equations. Finally, Section 5 presents 
the results of the numerical experiment, and Section 6 concludes the paper. 

2 General Rigid-Body Problem Statement 

For a generic rigid body the Euler's rotational equations of motion are (Gold- 
stein 1980) 

IiP = {h - h)qr + mi 

hq = {h - h) r p + m2 (1) 
hi" = (-^1 - h)pq + m3, 

where the dot symbol denotes the time derivative, {li : i = 1,2,3} are the 
moments of inertia about the three orthogonal body-fixed principal axes of 
inertia, {p, q, r} are the components along those axes of the angular velocity 
of the body with respect to the inertial frame and {nii : z = 1, 2, 3} are the 
components along the same axes of the external torque acting on the body. 

The rotation transformation matrix Rnb G 50(3) from a body fixed 
coordinate system B to the inertial coordinate system is defined such that 

ajv = Rnb o-b, (2) 
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where is the column matrix containing the components of a vector a in the 
Cartesian coordinate system A^, and as is the column matrix of components 
of the same vector in the Cartesian coordinate system B. 

The matrix Rnb obeys the following differential equation (Leimanis 
1965) 

Rnb = Rnb^, (3) 

where 

-r q 

r -p , (4) 

-q p 

Eq. 3 is the kinematic differential equation for the rigid body. Integration of 
this equation provides the current body orientation, and, together with the 
integrals of Eq. 1, completely solves the rigid body motion problem. 
In general. Equations 1 and 3 do not have exact analytic solutions. 



2.1 Kinematic differential equations in terms of complex stereographic 
variables 

Besides the rotation transformation matrix, several other set of variables can 
be used to describe the orientation, or attitude, of a rigid body (Shuster 
1993). 

In particular, the set of three complex stereographic rotation variables 
proposed in Longuski and Tsiotras 1995 is used in the present paper. The 
definition of these variables is summarized here below. 

Because Rnb is orthonormal, the elements of its k-th row (r^i, rk2 and 
Tfcs) are the coordinates of a point on the unit sphere in i?"^. 

The stereographic projection from the south pole of the unit sphere onto 
the equatorial complex plane establishes a map between and the extended 
complex plane CIJ oo. 
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In particular, the three points {(rfei, rj.2, ^'Ajs) : k ~ 1,2,3} are projected 
onto the following three complex numbers (Longuski and Tsiotras 1995) 

u>fc = — — , A: = l,2,3, (5) 

1 + 

where i denotes the imaginary unit. These constitute the stereographic rota- 
tion variables. 

Given the stereographic rotation variables, the elements of the matrix 
Rnb can be found by exploiting the inverse stereographic map from the 
complex plane onto the unit sphere. In particular, it results (Longuski and 
Tsiotras 1995) 

l + \wkr l + |u'fc|'^ l + \wk\^ 

where | • | denotes the absolute value of a complex number and the over bar 
indicates the complex conjugate. 

By differentiating Eq. 5 and by using Eq. 3, after some algebraic steps, 
the rotational kinematic differential equations in terms of the stereographic 
rotation variables arc found to be 

Wk^^ {p~iq) wl~irwk + ^ (p + iq), fc = l,2,3. (7) 

These equations, which are of the Riccati form, do not have in general an 
exact analytic solution. 



3 Case of a rigid body with spherical ellipsoid of inertia and 
constant torque 

For a rigid body having a spherical ellipsoid of inertia, i.e. either a uniform 
sphere or any uniform regular simplex (Klamkin, 1993), 



(8) 
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In this case, Eq. 1 become 

p = m\ 

q = m'^ (9) 
r = TO3, 

where {m[ = rrii/I : i = 1, 2, 3}. 

Let us now assume that the torque acting on the body is constant with 
respect to a body-fixed reference frame. 

A suitable rotational transformation of coordinates can be used to align 
the third axis of a new body-fixed coordinate system along the direction of 
the acting torque. Therefore, without loss of generality, Eq. 9 can be written 
as 

p = 

q = (10) 

r = [/, 

where p, q and r are used, from this point on, to indicate the components 
of the angular velocity along the new reference axes, and U indicates the 
constant torque, normalized by the value of the moment of inertia. 
The integration of Eq. 10 immediately gives 

p{t) = Po 

q{t) = qo (11) 
r{t) = ro + Ut, 

where po,qo and denote the initial condition values, and t denotes the 
time. 

Even if the solution of the dynamic problem has the simple linear form of 
Eq. 11, the corresponding kinematic differential equation, which is obtained 
by substituting the relations in Eq. 11 into Eq. 4 and 3, cannot be directly 



8 



integrated, with the exception of when the initial angular velocity is along 
the same axis of the acting torque {pq ~ qo = 0). In this case the solution is 

Rnb (t) = Rnb (0) exp (^J^ 12(0 , (12) 

where exp(») indicates the matrix exponential. Eq. 12 is only valid if the ma- 
trix f2{t) commutes with its integral (Adrianova 1995). This commutativity 
is lost when the initial angular velocity is arbitrary; therefore, Eq. 12 is not 
valid anymore in that case. The limitation in the validity of the solution of 
Eq. 12 is not clearly stated in Leimanis 1965, as already noticed in Longuski 
1984. 

In order to find the analytic solution for an arbitrary initial angular ve- 
locity, we consider the kinematic differential equations in terms of the stere- 
ographic variables. In particular, by substituting the relations in Eq. 11 into 
Eq. 7, it results 

w{t) = ^ (Po - iqo) w{tf -i{rQ + Ut) w{t) + ^ (po + no) (13) 
where the subindex k has been dropped in order to unburden the notation. 



4 Exact analytic integration of the kinematic differential equation 

Theorem 1 Given pQ,qQ,rQ and U real numbers, with Pq and Qq not both 
zero, the general solution of Eq. 13, governing the rotational kinematics of a 
rigid body with spherical ellipsoid of inertia, initial angular velocity compo- 
nents p{0) ~ po , q{0) = qo, and r{0) = ro along the three body fixed axes (Ti,(T2 
and (73, and subjected to a constant torque U, normalized by the value of the 
moment of inertia and directed along the axis a^, is the following 



with 



«^W = |^^^[6^ + G(z)], (14) 
•J [Po — iqo) 



(15) 



^ 2 iFi (g^, I; z^) (^ _ 1) ^2 + 6 iFi (1 - I, f ; z^) c^z - 3 (i^, f ; z^) 
iFi(-|,i;z2)c+iFi(if^,|;z2), 
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where iFi denotes the confluent hypergeometric function (Lebedev 1965), c € 
C is the constant of integration and 

{l + i){ro + Ut) ^ i{pl + qi) ^jg^ 



2VU 4[/ 
Proof Three successive variable substitutions are used in this proof. In the 
first substitution, the dependent variable in Eq. 13 is changed from w(t) to 
x{t) by the Riccati transformation (Zwillinger 1989) 

w{t) = - , ^' . . , 17 

x[t) [Po - m) 

then the equivalent second order homogeneous linear differential equation is 
obtained 

x{t) + i (ro + Ut) x{t) + sox{t) = 0, (18) 

where sq := {pi + ql) /4. 

In the second substitution, the dependent variable in Eq. 18 is changed 
from x{t) to y{t) by 

x(i) = y(t)exp -it(^Q + ^^ , (19) 

then the following equivalent equation is obtained 

y{t) - I (ro + Ut) y{t) + (sq - lU) y{t) = 0. (20) 

Finally, the independent and dependent variables in Eq. 20 are changed, 
from t to z, as defined in Eq.l6, and from y{t) to 77(2) = y{t) respectively. 
Therefore, the following equivalent equation is obtained 
d^riiz) „ dri(z) „ , ^ „ 

^-2z^+2.,(.).0, (21) 

where v is defined as in Eq. 16. This is a second order differential equation 
of Hermite type (Zwillinger 1989) and has the following general solution 
(Lebedev 1965) 

viz) = iFi (^-^, i; z^) a + ^F^ (^1^, | z^) bz, (22) 
with a e C and 6 G C constants of integration. 
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By substituting the expression in Eq. 22 baek into Equations 19 and 17, 
by defining c ~ a/b and by rearranging, the result stated in Eq. 14 and 16 is 
finally obtained. 

Corollary 1 An alternative expression of the solution given in Eq. 14 and 16 

is 

where H denotes the Hermite function (Lehedev 1965), c' ^ C is the constant 
of integration, and z and v are defined as in Eq.l6. 

Proof Since the Hermite function is by definition a particular solution of the 
Hermite equation, we can write the general solution of Eq. 21 as 

riiz) = iFi (^-^, ^; z^^ a' + H, (z) 6', (24) 

with a' e C and &' G C constants of integration. 

By using the expression in Eq. 24 instead of the one in Eq. 22, by defining 
c' = a' /b' , and by proceeding as for the proof of Theorem 1, the result in 
Eq. 23 is obtained. 

The confiucnt hypergeometric function used in Theorem 1. which is de- 
fined in C^, is (Lebedev 1965) 

OO - \ 

iFi(a,/3;7)=^i^^, H < oo, /? ^ 0, -1, -2, ... , (25) 
where (•)fc denotes the Pochammer symbol 

(•)o-l' (•)fc-fc(. + l)...(. + fc-l). A; = 1,2,.... (26) 



The series in Eq. 25 converges for all finite 7, and it represents an entire 
function of 7, for fixed a and (3. 
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4.1 Solution in terms of rotation matrix 



The problem of computing the rotation matrix function R^sit) is solved 
by finding the three complex functions wi{t),W2{t) and w^it), particular 
integrals of Eq. 13 which substituted into Eq. 6 give the elements of the first, 
second and third row, respectively, of the rotation matrix. 

Without loss of generality, we can consider the rotation matrix R^b 
be equal to the identity matrix at the initial time t = 0. Therefore, from 
Eq. 6, the initial conditions for wi{t),W2{t) and 'W3{t) are found to be 



(r-ii = l,ri2 = 0,ri3 = 0) ^ ^1(0) -i 
(r2i = 0,r22 = l,r'23 = 0) ^ ^2(0) = 1 
(^31 = 0, r32 = 0, r33 = 1) ^ W3(0) = 0. 



(27) 



The values {ci, 02,03} of the constants of integration, which substituted in 
the general integral (Eq. 14) give the three particular integrals wi (t) , W2 (t) 
and W3 (t) , are found by solving Eq. 14 for c and by imposing the three initial 
conditions of Eq. 27 



Cl 



C2 



C3 



(1 + *) 



6VU 



2 ' 2 ' 2f7 
V 3. i'': 



iFi ( i^, §; ^ (6r2 + 3 (po - «go) ''o + &iU) + 2g 



2 iFi (1 - f , i; ^ j i^ro + iFi (-f , i; g) {po - iqo + 2ro) 
iFi {6irl-3ipo~iqo)ro-6U)+2gi 



6VU 
(1 + *) 



I) i^ro + iFi (-f , |; i'Po + 10 + 2ro) 



u 3 . '"'1 
2 ' 2 ' 2(7 



2iFi [1^ 

.9 + 3iFi 



l-u 3. ^ 
2 ' 2' 2!7 



{rl + ^U) 



2^2^2Uj ^^^\ 2'2'2(7 

where 1^ is defined as in Eq. 16, and 

3 — 1/ 5 ivQ 



(28) 



iFi 



(^-l)^o- 



(29) 



2 '2'2[/^ 

An alternative method to compute the rotation matrix consists of using 
the above procedure to compute the elements of two rows of the matrix 
only, and then exploit the orthonormality conditions in order to find the 
elements of the remaining row. In particular, this method allows to overcome 
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the problem of the point at infinite of the stcreographie projection and the 
related singularity of the inverse stereographic map. In fact, whenever one of 
the body axes approaches the direction of the south pole of the unit sphere 
and, consequently, the correspondent point, mapped by the stereographic 
projection (see Eq. 5), tends to infinite, it is sufficient to consider only the 
remaining two axes in order to build the rotation matrix. This procedure 
makes the proposed analytic solution valid for any length of time and rotation 
amplitude. 

5 Numerical experiments 

A numerical experiment has been performed by using the newly found so- 
lution and its results are reported here below. The experiment considered a 
rigid body having spherical ellipsoid of inertia and moving with the following 
sample initial conditions and body fixed torque value 

Pa = 10, qa = 15, tq = 20, U ^ 3. (30) 

The numerical evaluation of the analytic solution is compared to the 
solution obtained by a numerical integration method at a sample particular 
time of t = 40. 

The numerical evaluation of the exact analytic solution in Eq. 14, with 
the constants of integration evaluated as in Eq. 28, gives the following values 
for the stereographic rotation variables 

wi = -0.4263512625196356 + 0.4033450082324547 i 

W2 = -0.2092049992603608 - 0.5057428050364536 i (31) 

W3 = 0.4167055594784583 + 0.1095946611492525*. 

By substituting the above values into the relations in Eq. 6, the corresponding 
rotation matrix results 



(32) 
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Rnb = 



-0.6000092673712773 -0.6342329852754623 0.4875832231087923 
0.7783397597095152 -0.3219671485837583 0.5390031295717849 
-0.1848677838995137 0.7029122815980806 0.6868320222985118 



The software package Mathematica was used to conduct the above numerical 
evaluations. This computation took 0.8 seconds of CPU time. 

For vahdation purposes, a direct numerical integration of Eq. 3 was con- 
ducted by using the Matlab algorithm odell3 (Adams-Bashforth-Moulton), 
with relative and absohitc tolerances set to 10"^*. This computation took 
11 seconds of CPU time. The numerical evaluation of the analytic solution 
given in Eq. 32 and the result of the numerical integration coincides up to the 
11*'' digits. This relative error between the two computations is attributed 
to the accuracy of the numerical integrator being lower than the one of the 
numerical evaluation of the analytic solution. This is confirmed by the cor- 
responding different level of accuracy in maintaining the orthonormality of 
the rotation matrix. 



6 Conclusions 

For the first time, to the knowledge of the author, the exact analytic solution 
has been introduced in this paper for the rotational motion of a rigid body 
having spherical ellipsoid of inertia (e.g. an homogeneous cube or sphere, or 
a sphere made of homogeneous concentric layers) subjected to a constant 
torque and arbitrary initial angular velocity. Therefore, the present paper 
adds a new significant element to the set of few special cases for which an 
exact solution of the rotational motion of a rigid body exists. Furthermore, 
the proposed analytic solution is valid for any length of time and rotation 
amplitude. 

In particular, the parametrization of the rotation by three complex stere- 
ographic variables was used in order to obtain a kinematic differential equa- 
tion of Riccati type which was then solved analytically by substitutions of 
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variables. The solution involves confluent hypcrgeomctric functions. The ro- 
tation matrix can be recovered from the complex rotation variables by inverse 
stercographic map. 

The results of the numerical experiment clearly confirms the exactness of 
the proposed analytic solution. 
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